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Abstract 

We searched for evidence of a Higgsstrahlung process in a secluded sector, leading 
to a final state with a dark photon U and a dark Higgs boson hi, with the KLOE 
detector at DA<hNE. We investigated the case of hi lighter than U, with U decaying 
into a muon pair and hi producing a missing energy signature. We found no evidence 
of the process and set upper limits to its parameters in the range 2 < mjj < 
1000 MeV, rrih' < mu- 
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1 Introduction 


Astrophysical data reveal in a more and more convincing way that our knowl¬ 
edge of the Universe is limited to about 4-5% of the total matter-energy con¬ 
tent: this is generally interpreted as an evidence of the existence of dark matter 
and dark energy components. In recent years, several astrophysical observa¬ 
tions have failed to find a common interpretation in terms of standard as¬ 
trophysical or particle physics sources (m 3 . Although there are alternative 
explanations for some of these results, they could all be explained with the ex¬ 
istence of a dark matter weakly interacting massive particle, WIMP, belonging 
to a secluded gauge sector under which the Standard Model (SM) particles are 
uncharged (111-201]. In a minimal model, a new abelian U(l)s gauge held is in¬ 
troduced, the U boson or dark photon, with mass near the GeV scale, coupled 
to the SM only through its kinetic mixing with the SM hypercharge held. The 
kinetic mixing parameter e is expected to be of the order 10 -4 — _1CV 2 fl2 21 ]. 
so that observable effects can be detected at e + e _ colliders 121-251 or at hxed 


target experiments working in the GeV region [261-29]. The existence of the U 
boson, through its mixing with the ordinary photon, can also accomodate the 
observed discrepancy in the measured muon anomalous magnetic moment a^ 
with respect to the SM prediction jdoj. Several searches of the U boson have 
been performed in recent years with negative results, setting upper limits to 
e: Al [3d!, 0, APEX 0, WASA 
3 §. 


HADES 35], KLOE |36|, [37], BaBar 


Since the U boson needs to be massive, one can implement, in close analogy 
with the SM, a spontaneous breaking mechanism of the 17(1)5 symmetry, thus 
introducing a Higgs-like particle, b! or dark Higgs, whose mass hierarchy with 
the dark photon is not constrained by the theory [0 . 

The U boson can be produced at e + e _ colliders via different processes: e + e~ —> 
f/y , e + e _ —* Uh' (dark Higgsstrahlung) and in decays of vector particles to 
pseudoscalars. In this work the Higgsstrahlung process e + e _ —* Uh' is studied, 
using data collected by the KLOE experiment at the e + e _ collider DA<hNE 
at the Frascati laboratory, both at a center of mass energy of ~ 1019 MeV, 
the mass of the 0 meson (on-peak sample), and at a center of mass energy of 
~ 1000 MeV (off-peak sample). The process e + e _ —>■ Uh', with U decaying 
into lepton or hadron pairs, is an interesting reaction to be studied at an e + e~ 
collider, being less suppressed, in terms of the mixing parameter, than the 
other final states listed above. There are two very different scenarios depending 
on the masses of the dark photon mjj and of the dark Higgs boson m//. For m^/ 
larger than 2 mu, the dark Higgs boson would decay dominantly and promptly 
to a U boson pair, thus giving rise to a six charged particle final state (the 
scenario with mw larger than mu but smaller than 2 mu is similar, with one 
dark photon off shell): this case was recently investigated by the BaBar [39] 
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and Belle (4o| experiments. On the other side, Higgs bosons lighter than the 
dark photon would have, in most of the parameter space region, such a large 
lifetime to escape detection, showing up as a missing energy signature. We 
confined the search only to the latter case, m/,/ < mu, the so called “invisible” 
dark Higgs scenario. 


The lifetime of the dark Higgs boson depends on the kinetic mixing parameter 


e, the boson masses m^ and mu and the dark coupling constant a.£> 22], For 
boson masses of the order of 100 MeV and an = oi ern , the dark Higgs boson 
lifetime would be ~ 5/is for e ~ 10 -3 , corresponding, for the energy range 
explored in this analysis, to a decay length of ~ 100 m. The dark Higgs boson 
would be thus invisible up to e ~ 10~ 2 -=- 10 -1 , depending on the h' mass. 


In this work the search is limited to the decay of the U boson in a muon pair: 
the final state signature is then a pair of opposite charge muons plus missing 
energy. The measurement is thus performed in the range 2m fl < mu < 1000 
MeV with the constraint my < mu- 


The production cross section of the dark Higgsstrahlung process is propor¬ 
tional to the product a^xe 2 and depends on the boson masses 22], Values as 


high as hundreds of fb are reachable in this model. Compared to the B-factory 
case j3Si, 40], KLOE benefits of the 1/s factor and of the resonance-like be¬ 
haviour expected for the production cross section [ 22 J . The branching ratio of 
the U boson into muon pairs is predicted to be just below the 50% level for 
masses slightly above the kinematical threshold mu = 2m M , then to decrease 
up to a minimum around 5%, for masses corresponding to the p resonance (due 
to the concurrent decay into hadrons), and then to increase to ~ 30 -=-40% up 
to mu — 1 GeV 221 . 


2 The KLOE Detector 


DA$NE, the Frascati 4>-factory, is an e + e _ collider working at the center of 
mass energy, y/s ~ m§ = 1.0195 GeV [41]. Positron and electron beams col¬ 
lide at an angle of 7r—25 mrad, producing c{) mesons nearly at rest. The KLOE 
detector is made up of a large cylindrical drift chamber (DC) [421. surrounded 
by a lead scintillating fiber electromagnetic calorimeter (EMC) [43[ . A super¬ 
conducting coil around the EMC provides a 0.52 T magnetic held along the 
axis of the colliding beams. 

The EMC consists of barrel and end-cap modules covering 98% of the solid 
angle. The calorimeter modules are segmented into five layers in depth and 
read out at both ends by 4880 photomultipliers. Energy and time resolutions 
are cte/E = 0.057/ yjE(GeV) and a t = 57ps /yjE{ GeV) © lOOps, respec- 
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tively. The drift chamber, with only stereo sense wires, 4 m in diameter and 
3.3 m long, has a mechanical structure in carbon-fiber and operates with a 
low-mass gas mixture (90% helium, 10% isobutane). The spatial resolutions 
150 p.m and a z ~ 2 mm. The momentum resolution for large angle 


are a X y rs -' 


tracks is a p± /p± ~ 0.4%. The trigger 44| uses both EMC and DC information. 


Data are then analysed by an event classification filter 45| , which selects and 


streams various categories of events in different output hies. 


3 Event selection 


The analysis of the process e + e~ — >■ Uh' , U —> /i + /r~, hi invisible (e + e _ —> 
Uh' in the following), has been performed on a data sample of 1.65 fb _1 col¬ 
lected at a center of mass energy of ~ 1019 MeV, corresponding to the mass 
of the 0 meson (on-peak sample in the following), and on a data sample of 
0.206 fb _1 at a center of mass energy of ~ 1000 MeV (off-peak sample in the 
following), well below the 0 resonance. 



COS0 


Fig. 1. Distribution of the polar angle of the muon pair momentum for the signal 
e + e _ —>• Uh' (black line) and for e + e~ —>• 7 (red line). Here the two processes 

are not normalised and are shown only in order to compare the shapes of the dis¬ 
tributions. All the generated samples at various mu and mw are included in the 
signal distribution. 

The Monte Carlo simulation of the signal process e + e _ —y Uh' has been pro¬ 
duced using an ad hoc generator interfaced with the standard KLOE simula- 
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tion program |45|. Signal samples have been generated for various pairs of mh> 
- mu values along a grid with steps of ~ 30 MeV to cover all the allowed kine¬ 
matic region. The invariant mass resolution varies between 0.5 and 2 MeV for 
the muon pair (M^), and between 3 and 17 MeV for the event missing mass 
(M miss ). The signal process signature would thus be the appearance of a sharp 
peak in the bidimensional distribution - M miss . Moreover, the distribu¬ 
tion of the polar angle direction of the muon pair momentum, 9 , contrarily 
to most of the dominant background processes, is expected to prefer large 
angles. The differential cross section has two dominant terms proportional to 
sin 9 and sin 3 9 j22|, with relative weights smoothly dependent on the boson 
masses. This angular distribution allows to reject most of the background of 
QED processes with a simple geometrical selection and implies that the miss¬ 
ing momentum direction preferably points to a very well equipped region of 
the KLOE detector, where the best efficiency is achieved. Fig.|T] shows the 
distributions of the muon pair polar angle direction for the signal e + e^ —y Uh' 
(black line) and the e + e _ —> /i + /r _ 7 background (red line), where all the 
generated samples at various mu and rriu masses are included in the signal 
sample. 

As a first step of the analysis, a preselection was performed by requiring: 

• events with only two opposite charge tracks with associated EMC clusters, 
with polar angles | cosffi^l < 0.8 and momenta below 460 MeV, that form a 
reconstructed vertex inside a cylinder of 30 cm length, 4 cm radius, centered 
at the interaction point (IP); 

• the sum of the momenta of the two tracks to be greater than 450 MeV; 

• the polar angle of the dimuon momentum to be in the barrel acceptance: 

| cos^l < 0.75; 

• the modulus of the missing momentum to exceed 40 MeV. 

After this selection, mostly aimed at rejecting backgrounds from QED pro¬ 
cesses, the hermeticity and tightness of the electromagnetic calorimeter was 
used as a veto to avoid the presence of photons in the event by requiring 
no prompt EMC clusters unassociated to tracks. The calorimeter veto inef¬ 
ficiency as a function of the energy was studied with a sample of radiative 
Bhabha scattering events e + e _ —> e + e _ 7 and found to range between 10% at 
20 Mev and 0.1% at about 200 MeV. 

The event selection then proceeded by applying a particle identification (PID) 
algorithm to the two tracks, based on the excellent energy and time resolution 
of the EMC. A set of feed-forward neural networks, organised for different 
values of track momentum and track polar angle, was trained on simulated 
Monte Carlo samples to perform muon to electron discrimination. The neural 
networks used five input variables (cluster time, energy to momentum ratio 
and three variables related to energy depositions in calorimeter layers) and 
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produces one output. The PID performances were checked on selected data 
samples of e + e _ —* e + e _ , e + e^ —> e + e _ —* 7r + 7r _ : the fraction of 

events where both tracks were identified as muons was measured to be 85% 
in e + e _ —* /r + /i - events, 10 -4 in e + e _ —* e + e _ events and ~ 50% in e + e _ —* 
7r + 7r _ events (showers produced by muons or pions have similar properties at 
low energies). 

After missing energy and PID selections, a large background from 0 —>■ K + K~, 
K ± — > fi^u events survived in the on-peak sample. This happens when both 
kaons decay semileptonically close to the IP. Charged kaons have an average 
decay length of ~ 90 cm in KLOE. The reconstructed vertex of the muon 
tracks is thus expected to be displaced from the IP and with a bad fit quality. 
Cuts on the radial and z projections of the distance between the reconstructed 
vertex and the IP and on the y 2 of the fit allowed to reduce by a factor ~ 80 
the 0 —> K + K~, K± —> background, lowering the signal efficiency by ~ 

65%. 

Events surviving all the described selections were organized in bidimensional 
histograms with the muon pair mass and the event missing mass M miss 
on the two axes. The binning was chosen to keep most of the signal inside a 
single bin. For a 5 MeV bin width was enough over all the plane; while for 
Mmiss a variable binning of 15, 30 and 50 MeV widths was chosen. According 
to the simulation, a fraction of 90-y95% of the signal was contained in one 
single bin. The signature of the process would thus be the appearance of an 
excess in a single bin in the M^^-M m i SS plane over the background. The signal 
selection efficiency, estimated from Monte Carlo on the generated points of 
the mu-mh' grid, was found to be between 15% and 25%, depending on the 
masses, with most frequent values of ~20%. The efficiency for a generic point 
on the M^^-Mmiss plane was then evaluated by linear interpolation. 


4 Results 


After all the described selections, 15278 events survived in the on-peak sample 
(fig. [21 left plot) and 783 in the off-peak sample (ffg.EJ right plot). In the left 
plot of fig. [2] (on-peak sample) several sources of backgrounds can be distin¬ 
guished: 

• 0 —> K + K~, K ± —y ^v (quadrangular region at the left of the populated 
part of the distribution); 

• (j) —> n + n~n° (quasi-horizontal band, corresponding to events in which both 

photons from the 7r° decay are undetected), partly intersecting the 0 —* 
Ii + K~, K ± —>■ region; 

• e + e _ —>■ /r + /i _ and e + e _ —> 7r + n~ events in the continuum (diagonal and 
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horizontal bands starting from the right-bottom part of the distribution); 

• e + e _ —* e + e - /r + /i~ and e + e~ —>■ e + e _ 7 r + 7 r^ (photon-photon interactions, 
top triangular part of the distribution, for M miss > 350 MeV), with e ± in 
the final state being scattered at very small angles in the beam pipe. 

In the distribution in the right plot of fig. [2] (off-peak sample) all the back¬ 
grounds from the </> decays are strongly suppressed and only those in the 
continuum remain visible. 
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Fig. 2. Results for on-peak sample (left plot, 1.65 fb 1 integrated luminosity) and 
off-peak sample (right plot, 0.206 fb -1 integrated luminosity). 


Monte Carlo generators fully interfaced with the KLOE detector simulation 
program were available for all the background processes but for the e + e _ —y 
e + e~fi + n~ and e + e _ —> e + e~7i + 7T~. For these two processes the Courau gen¬ 
erator program 46J] was used and the results smeared to keep into account the 
detector effects. 


As most of the signal is expected to populate a single bin of the mass dis¬ 
tributions, a 5x5 bin matrix in the M^-M miss plane was built and moved 
sliding all along the distributions of fig. [2] both on data and Monte Carlo. The 
presence of a possible signal was checked by using the central bin, while the 
others were used for background evaluation. This was done by computing a 
data-Monte Carlo scale factor k based on the sum of the contents of the 24 
bins surrounding the central one in data {DT 2 a) and Monte Carlo (MC 24 ): 
k = YTctl' The prediction for the background in the central bin is then simply 
defined as the product of the central bin content in Monte Carlo rescaled by k. 
The usage of the described scaling procedure allowed to reduce the systematic 
uncertainties due to the background evaluation (see sec. 5). Fig. [3] shows the 
data-Monte Carlo comparison after the scaling correction for the on-peak and 
off-peak samples, projected along the M /t/i and M miss axes, together with the 
individual contributions of the different background processes. The agreement 












Fig. 3. Data - Monte Carlo comparison for the on-peak sample (top plots) and 
off-peak sample (bottom plots). Projections along the axis (left plots); projec¬ 
tions along the M m i ss axis (right plots). Also shown are the various contributing 
backgrounds. 

is satisfactory all over the populated regions of the distributions. 


5 Systematic errors 


Systematic uncertainties affect the signal efficiency evaluation and the back¬ 
ground estimate. Several sources of systematic uncertainties in the signal effi¬ 
ciency evaluation from Monte Carlo were taken into account. 

Uncertainties from the PID procedure were estimated by selecting samples of 
e + e _ —* /r + /i “7 in data and Monte Carlo, applying the PID algorithms to 
one of the two tracks to increase the purity of the selection and studying on 
the opposite track the data-Monte Carlo differences of the PID efficiency as 
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a function of the track momentum. The total effect, defined as the average 
product of individual effects on the single tracks, was found to vary between 
2% and 3%, depending on the boson masses. A similar procedure was applied 
to evaluate the correction factors and systematic uncertainties of the PID 
algorithms for pion identification, which affect the background evaluation. 

The same e + e _ —» yU + /i _ 7 samples selected in data and in the simulation were 
used to evaluate the effect of the cut on the vertex-IP distance. A correction 
to the Monte Carlo signal efficiencies of the order of 15%, weakly dependent 
on cos 9, was derived and applied. An associated systematic error of 0.5% was 
estimated and added on the signal efficiency evaluation. 

The systematic uncertainty due to the usage of the EMC veto was evaluated 
by selecting samples of (f) —> K + K~, K ± —> in data and Monte Carlo. 

In this case, the cut on the vertex-IP distance was slightly relaxed, in order 
to increase the size of the sample. A 2% data-Monte Carlo difference was 
observed and used both to correct the Monte Carlo efficiency and to quote a 
systematic uncertainty due to this source. 

The systematic uncertainty due to the kinematical preselections of the analysis 
was estimated by varying track angles and momenta within their measurement 
errors by one standard deviation: a 1% effect was ascribed to this source. 

The systematic uncertainty due to the binning choice was estimated by eval¬ 
uating in the simulation the binomial statistical error on the fraction of the 
signal contained in one bin. This turned out to be of the order of 0.3%, on 
average. 

Finally, an average ~ 1% uncertainty was estimated due to the linear inter¬ 
polation procedure in the signal efficiency evaluation process. 

The total systematic uncertainty on the signal efficiency was then evaluated as 
the quadratic sum of all the above effects. It never exceeded 4%, with an av¬ 
erage value of 3.5%, very small when compared to the statistical uncertainties 
affecting this measurement. 

Most of the systematic uncertainties in the background evaluation cancel in 
the scale factor ratio k. All the systematic sources considered for the signal 
efficiency evaluation, but those related to the linear interpolation procedure, 
were taken into account and their effect on the background estimate computed. 

Additional effects were taken into account. The uncertainties on the back¬ 
ground process cross sections were varied within their theoretical and mea¬ 
surement errors; a further 1% uncertainty was added for those related to the 
photon-photon final states, for which no full simulation was available; the 
uncertainty on the integrated luminosity was estimated to be 0.3%. 
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The total systematic uncertainty on the background was evaluated as the 
quadratic sum of all the above effects. It has an average value of 2.7% with a 
very small tail extending up to 10%. 


6 p 0 values and Upper Limits 


In order to evaluate the compatibility of the observed results with the back¬ 
ground only hypothesis (p 0 value) and to derive upper limits to the parameters 
of the dark Higgstrahlung process, a Bayesian procedure was set up. For each 
position of the 5x5 bin mass matrix, a likelihood function was devised based 
on uniform prior probabilities of the counting variables (constrained to be non 
negative) and on four Poissonian distributions representing the probabilities 
related respectively to the number of observed events in the central bin of 
the sliding matrix, the number of predicted background events in the same 
bin from Monte Carlo, the number of observed and predicted events in the 
surronding 24 bins (DT 2 4 and MC 2 4 , entering in the scale factor ratio k). This 
procedure takes thus into full account the fluctuations due to the data and 
Monte Carlo statistics. The systematic uncertainties on the signal efficiency 
and on the background estimate were taken into account by convolving the 
four Poissonian distributions with two correlated Gaussian distributions, with 
variances set equal to the estimated systematic errors. Whenever the dark 
Higgsstrahlung process was searched for, the small fraction of signal expected 
outside the central bin of the 5x5 sliding matrix was explicitly taken into 
account in the likelihood expression. 



Fig. 4. Left: po value distribution for the on-peak sample (red line), off-peak sample 
(blue line), combined sample (black line). Right: po values for the combined result. 

The p 0 distributions for the on-peak, off-peak and combined samples are shown 
in fig.H left plot. Fig. SI right plot, shows the computed po values as a func- 
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tion of M^-Mmiss masses for the combined sample. There are three values 
exceeding the threshold corresponding to a 3a excess, while 4.2 were expected 
on probabilistic base. The excess significance of those points (see fig. HI) are 
at the level of 3.1<r, 3.2a and 3.4a. In the on-peak and off-peak samples the 
most significant values exceeding the 3a threshold are at the level of 3.9a and 
3.8a respectively (see fig. HI left plot) These excesses, even though at quite in¬ 
teresting level, are then lost in the combination of the two samples, becoming 
fluctuations of average size. 


As no evidence of the dark Higgsstrahlung process was found, 90% confidence 
level Bayesian upper limits on the number of events were derived bin by bin 
in the M IHt -M rniss plane, separately for the on-peak and off-peak samples, and 
then converted in terms of a# x e 2 . They are shown in fig. [5] Fig. [6] shows the 
on-peak and off-peak 90% CL upper limits projected along the m v and my 
axes after a slight smoothing to make them more readable. The different curves 
in mu (my) correspond to different values of my (mu)- These results were 
then combined by taking into account the different integrated luminosities of 
the two samples and the respective signal efficiencies and cross sections. The 
combined results are almost everywhere dominated by the on-peak sample, 
because of the larger available statistics, with the exception of some very 
noisy background regions. They are shown in fig. [7] These limits are largely 
dominated by the data statistics. Values as low as lCT 9 -ylCV 8 of the product 
«d x e 2 are excluded at 90% CL for a large range of the dark photon and dark 
Higgs masses. 


a D s 2 



100 200 300 400 500 600 700 800 900 1000 
my (MeV) 


a D s 2 



200 400 600 800 1000 

my (MeV) 


Fig. 5. 90% CL upper limits in au x e 2 for the on-peak sample (left plot) and 
off-peak sample (right plot). 
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Fig. 6. Top plots: 90% CL upper limits in old x e 2 as a function of mu for different 
values of mh 1 for the on-peak sample (top, left) and off-peak sample (top, right). 
Bottom plots: same limits as a function of my for different values of mu for the 
on-peak sample (bottom, left) and off-peak sample (bottom, right). 


7 Conclusions 


A search for the dark Higgsstrahlung process e + e _ —» Uh ', U —> /i + /r _ , h! 
invisible, has been performed by KLOE in the range 2 < mu < 1000 MeV 
with mh' < mu- No evidence for signal has been observed and upper limits on 
the product of the kinetic mixing parameter e and the dark coupling constant 
old have been set in the range 10~ 9 -bl0~ 8 in old x e 2 - With the arbitrary 
hypothesis old = OL em these measurements translate into limits on the kinetic 
mixing parameter e in the range 10~ 4 -bl0” 3 . 
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Fig. 7. Combined 90% CL upper limits in x e 2 as a function of mjj for different 
values of my (top plot) and as a function of my for different values of mu (bottom 
plot). 
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